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Abstract 

We derive the long- wavelength non-linear sigma model for a two-dimensional Heisenberg system in the presence of 
the Dzyaloshinskii-Moriya and pseudodipolar interactions. We show that the system is a non-conventional easy-axis 
antiferromagnet, displaying an anomalous coupling between the magnetic field and the staggered order parameter. 
Our results are in good agreement with recent experimental data for undoped La2Cu04 compounds. 
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It is well established that spin fluctuations play a 
prominent role in high-temperature superconductors. 
In particular, the parent compounds display an anti- 
ferromagnetic (AF) phase which is usually described 
as a two-dimensional (2D) Heisenberg antiferromag- 
net [1]. However, recent magnetic-susceptibility [2] and 
Raman [3] experiments in detwinned La2Cu04 (LCO) 
single crystals can not be understood within the theory 
of a conventional antiferromagnet. In Ref. [4] we ar- 
gued that these anomalies originate from the presence, 
in the low-temperature orthorhombic (LTO) phase of 
LCO, of the anisotropic Dzyaloshinskii-Moriya (DM) 
and pseudodipolar (XY) interactions between spins. 
In this paper we derive explicitly the long-wavelength 
limit of the following 2D S = 1/2 Hamiltonian: 

H = J2 JS < ' S ^ + D ^ ' ( Sl x S J') + S * r y' S J> (!) 

(id) 

where Si represent the Cu ++ spins, J is the AF super- 
exchange, and Dij and Tij are, respectively, the DM 
vectors and XY matrices [5]. In the (xyz) coordinate 
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system of Fig. 1 the DM vectors are oscillating be- 
tween neighboring bonds with Dab = (0,d, 0), and 
Dac = (d, 0, 0), while the F matrices are constant and 
given by Tab = diag(r x + r 2 , Ii - T 2 , T 3 ) and Y B c = 
diag(Fi - r 2 ,ri + r 2 ,r 3 ). Here d and Fi, 2 ,3 > are 
of order 10 -2 and 10 -4 , respectively, in units of J [5]. 
To construct the long wavelength effective theory for 
the above Hamiltonian we decompose the unit vector 
fii = Si/S at site ri into its slowly- varying staggered 
and uniform components, Q,i — e lC ^ ri rii + aLi, where 
Q = (ir, 7t) and a is the lattice parameter. The con- 
straint 0% = 1 is enforced by n 2 = 1 and Li • n; = 0. 
Using this decomposition, the Heisenberg part of the 
Hamiltonian (1) has the standard form [6] 




while the DM and XY terms are given by 



Com = Dij • (Si x Sj) 

= 1^! J d 2 x d+ • (n x L) (3) 

and 



Preprint submitted to Elsevier Science 



6 February 2008 



Cxy = SiTijSj = ~ [ d 2 x \2Tfn a n p + 
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a?T°f B (d x n a d x np) - \o 2 Y°^ c {d y n a dyn f3 ) 



(4) 



Observe that only the combination d+ = (Das + 
Dac)/2 = {d, d, 0)/2, oriented along the a direction 
in the LTO coordinate system, and T+ = (Fas + 
Fac)/2 = diag(ri, Ti, T$) enter the continuum limit 
of the model (1). Moreover, since r, < J we can ne- 
glect the anisotropic corrections in the second line of 
Eq. (4) to the (Vn) 2 term of Eq. (2). Using the fact 
that n 2 = 1 we can rewrite the first term of Eq. (4) as 
Cxy = const + (2S 2 /a 2 ) J d 2 x(r x - T 3 )n 2 , which ex- 
plicitly indicates that the properties of the model (1) 
only depend on the difference Fi — T3 > 0. To derive 
the Euclidean action of the model (1) we introduce the 
path-integral coherent states representation of the spin 
states, which in addition to the previous contributions 
gives rise to the (dynamical) Wess-Zumino term [6] 



Cwz = — — / d 2 xL ■ (n x h) 



(5) 



so that the partition function Z = J DnS(n 2 — l)e~ s , 
with the action 5 = J &t[Ch.i + Cdm + Cxy + Cwz\- 
After integration of the L fluctuations we obtain in 
the LTO coordinate system the following modified 
non-linear a model (NLSM) (/3 = 1/T and f = 
/cfdr/d^): 



s = 



± j {(drn) 2 +c 2 (X7n) 2 + D 2 + n 2 a +r c n 2 } . (6) 



The bare coupling constant g and spin velocity c are 
given by gc = 8Ja 2 and c = 2\/2JSa, and we defined 
D+ = ^/2gcS 2 /Ja 2 d+ = (2V2Sd)e a (e„ is the unit 
vector along the a direction) and T c — (AgcS 2 /a 2 ) (Ti — 
r 3 ) = 32JS ,2 (Fi - r 3 ). Eq. (6) describes a 2D anti- 
ferromagnet with an easy axis along the b direction. 
Moreover, in the ordered AF phase the system displays 
also an uniform magnetization along the c direction. 
In fact, the saddle-point value for L, calculated from 
the Euclidean action S, is 



8i/ (nXh)+ 2^ (nXd+) 



(7) 



Thus, when (n) = aoet we have (L) = (l/8aSJ)(croef,X 
D+e a ) = (a D + /8aSJ)e c . 

We can extend the previous result in Eq. (6) for the 
case of a multi-layered quasi-2D system, with an ex- 
change coupling J± along the c direction, by writing 
5 -> (l/2ffc) J2 m I £™, where 



where n m is the staggered magnetization in the rath 
layer. We note that for a layered system the vector 
D+ = (— l) m D+e a alternates in sign between neigh- 
boring layers, giving rise also to the staggered order 
of the uniform L m components, in agreement with the 
experiments [2,3]. 
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Fig. 1. Left: tetragonal (xyz) and orthorhombic (bac) reference 
system for the LCO plane. The hatched/empty circles repre- 
sent the O ions tilted above/below the CUO2 plane. Right: 
arrangement of the staggered magnetization (small black ar- 
rows) and DM vectors (open arrows). 

Finally, in the presence of an uniform magnetic field 
the total Lagrangian C = Chj + Cdm + Cxy + Cwz 
acquires an additional term Cb = —(S/a) J d 2 xB ■ 
[L m — (L m ■ n m )n m ], where we made explicit the con- 
straint L m • n m = 0. This leads to an extra contribution 
(l/8aJ5)[B-(B-n ro )n TO ] to Eq. (7), which after inte- 
gration over L m has two effects on Eq. (6): (i) the mod- 
ification 9 r n m — > <9 T n m +iB x n m of the time deriva- 
tive [7]; (ii) an additional term (1/gc) J B- (D+ x n m ), 
for each layer. Observe that while the XY interaction 
only provides an easy-plane anisotropy for the stan- 
dard Heisenberg system, via the r c n 2 term of Eq. (6), 
the DM interaction gives rise to an unconventional 
coupling between the uniform magnetic field and the 
staggered magnetization. As discussed in detail in Ref. 
[4], this coupling is responsible for all the magnetic- 
susceptibility anomalies observed experimentally [2]. 



References 

[1] see M. A. Kastner et al, Rev. Mod. Phys. 70, 897 (1998) 
and references therein. 

[2] A. N. Lavrov et al, Phys. Rev. Lett. 87, 017007 (2001). 

[3] A. Gozar et al, Phys. Rev. Lett. 93, 027001 (2004). 

[4] M. B. Silva Neto, L. Benfa tto, V. Juricic, and C. Morais 
Smith, |cond-m at/0502588 to appear in Phys. Rev. B 
(2005). 

[5] L. Shekhtman et al, Phys. Rev. Lett. 69, 836 (1992): W. 
Koshibae et al, Phys. Rev. B 50, 3767 (1994). 



2 



I. Affleck, J. Phys. Cond. Matt. 1, 3047 (1989). 
S. Allen, D. Loss, Physica A 239, 47 (1997). 



